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Abstract. Pyrochlore dielectric materials such as (Bii 5Zno.5)(Nbi 5Zno.5)07 (BZN) 
have generated interest because they combine high dielectric constants with small 
dielectric loss tangents and yet are cubic at all temperatures. The recent low- 
temperature preparation and structural characterization of Bi2Ti207, which remains 
cubic down to 2 K, has provided a good model system for understanding the 
properties of Bi-based pyrochlores. In this contribution, the electronic structure 
of cubic Bi2Ti207 is visualized and compared with the electronic structure of 
the Aurivillius phase ferroelectric SrBi2Ta209 (SBT), which displays a ferroelectric 
distortion below 608 K associated with the tendency of lone pair active Bi''+ to move 
off-center. Such coherent off-centering distortions are frustrated on the pyrochlore 
lattice, and this prevents a ferroelectric-paraelectric phase transition in Bi2Ti207. 
Instead, Bi'^"'" ions in Bi2Ti207 are obliged to off-center in an incoherent manner, that is 
compatible with the cubic structure being retained. Frustrated lone pair behavior in 
the defect pyrochlore Pb2Sn206 is also described. Parallels between the well-studied 
frustration of certain types of magnetism in pyrochlore compounds (spin-ice) and 
the striking paucity of ferroelectric pyrochlores, arising from the corner-connected 
tetrahedral topology of the pyrochlore lattice are pointed out. 




Figure 1. (Color) (a) Ideal pyrochlore A2B2O6O' crystal structure in the Fci3m space 
group showing black A atoms, orange O' and the network of corner-connected BOg 
octahedra (blue), (b) is a view of the tetrahedral A2O' sublattice, and (c) is a view of 
the B2O6 sublattice. All views are looking down [110]. 



1. Introduction 

Pyrochlore oxides[l] have attracted a great deal of recent attention as a consequence 
of their diverse, often unusual, physical properties. Recent examples of the 
diverse properties displayed by pyrochlores include the discovery of colossal 
magnetoresistance in TI2M02O7 in the absence of multiple valence states,[2] 
superconductivity in Cd2Re207,[3] and in the /3-pyrochlore KOs206,[4] and an 
unusual chiral magnetic ground state in Nd2Mo207.[5] 

In the pyrochlore crystal structure, originally associated with the mineral 
Ca2Nb206F, the B (Nb) atoms form a network of comer-sharing octahedra. The A- 
site (Ca) has six O neighbors and two F neighbors. For the purpose of this work 
on oxide pyrochlores, the structure is best described as A2B2O6O' comprising two 
interpenetrating sub lattices with the formulae A2O' and BOg. These views of the 
structure are displayed in figure. 1. The crystal chemical principles involved in the 
pyrochlore structure are unusual as has recently been emphasized by Vanderah and 
coworkers. [6] For example, the A2O' sublattice can tolerate a great deal of displacive 
disorder, and consequently, does not obey the usual rules for cation substitution. 
These unusual principles play a significant role in rendering the pyrochlore structure 
suitable for the storage of ions involved in radioactive waste. [7] 

Another crucial feature of the structure is that when regarded separately, both 
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the A and B sublattices form comer-connected networks of A4 and B4 tetrahedra. The 
formation of such tetrahedral networks can result in magnetic frustration associated, 
for example, with the tnfeasibility of decorating the vertices of a tetrahedron with all 
four spins in an anti-parallel configuration. [8, 9] P. W. Anderson[10] first recognized 
the relation between the tetrahedral spin configurations of spinels and the problem of 
the crystal structure of ice-//,. According to the Bemal-Fowler rules[ll] each oxygen 
atom in the ice lattice is tetrahedrally bonded to four protons, with the proviso that 
two of the protons are proximal (covalent) and two are distal (H-bonded). There 
is no unique arrangement which satisfies these rules, and as a consequence, ice-7/, 
has a residual zero- temperature entropy [12, 13, 14] The crystal structure of ice- J/, is 
therefore frustrated in the sense that it cannot find a perfectly crystalline zero-entropy 
ground state. There is great current interest in such frustrated systems, associated 
with a desire to understand the nature of complex potential energy landscapes, 
and a number of recent studies have emerged on "spin-ice" pyrochlores.[15, 16] 
Interestingly, when ground state structures are induced, for example through OH~ 
doping, polar ground states become possible in water ice. [17] 

Unlike ABO3 perovskite oxides where ferroelectricity is not uncommon and is a 
well-studied phenomenon,[18] there exist very few ferroelectric pyrochlores. The 
few known system seem to be based on Cd2Nb207,[19] and its variants. [20, 21] 
These compounds remain the subject of study. [22] "High-fc" pyrochlore oxides with 
large dielectric constants exemplified by (Bii,5Zno.5)(Nbi.5Zno.5)07 (BZN) however 
do exist and are attracting considerable attention. BZN is a cubic pyrochlore 
at all temperatures,[23] with considerable short-range order in the structure. [24] 
Ceramic samples are known to combine high dielectric constants with low dielectric 
losses. [25, 26] A combination of spectroscopic techniques reveal that the unusual 
properties of BZN arise due to local hopping of the Bi and O' atoms among 
several potential minima. [27] Interestingly, from the viewpoint of applications, these 
desirable characteristics are retained in thin films. [28, 29, 30, 31] 

It is intriguing that BZN exhibits all the characteristics of a compound about to 
undergo a phase transition, such as the need to split the Bi and O' for satisfactory 
Rietveld refinements of neutron diffraction data,[23] and yet remains cubic till the 
lowest temperature. The absence of a phase transition is an important ingredient in 
the usefulness of BZN. For example, in a material with no phase transition, there is 
little fatigue, and the temperature coefficient of the dielectric constant is small. 

The preparation of a slightly off-stoichiometric pyrochlore phase related 
to Bi2Ti207 by Sleight and coworkers, [32] and the more recent preparation 
of stoichiometric, cubic Bi2Ti207 by Hector and Wiggin[33] provide a good 
model system for understanding polar phenomena in the bismuth pyrochlores. 
Stoichiometric Bi2Ti207 has been prepared by a low-temperature route and 
established by powder neutron diffraction to be a cubic pyrochlore down to 2 K, 
but with extensive site disorder at both the Bi and O' sites. [33] Site disorder on 
the pyrochlore A-site is well-known for lone pair active cations such as Tl+,[34] 
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Sn^+,[35, 36] and 51"^+ [23, 32, 37] and will be discussed in further detail. 

We have for some time been interested in the stereochemistry of lone pairs in 
extended solids,[38, 39, 40, 41, 42] and have used density functional calculations 
of electronic structure in conjunction with the electron localization function[49, 50] 
for visualizing lone pair electrons in the real space of crystal structures. Here such 
analysis is extended to the electronic structure of cubic Bi2Ti207. For comparison, 
LMTO calculations have been performed on the Aurivillius phase ferroelectric 
SrBi2Ta209[43, 44] where the transition from a paraelectric to a ferroelectric phase 
at 608 K is driven in large part by the off-centering tendency of the Bi'^+ ions. 
The electronic structures of paraelectric [45] and ferroelectric SrBi2Ta209[46, 47] have 
been described previously, but no discussion of the disposition of lone pairs in the 
ferroelectric phase has been presented. The disposition of Pb^^ lone pairs in the 
cubic defect pyrochlore Pb2Sn206[48] are also examined. 

2. Details of computation 

Linear Muffin-Tin Orbital calculations[51, 52] on pyrochlore Bi2Ti207 were 
performed within the atomic sphere approximation using version 47C of the Stuttgart 
TB-LMTO-ASA program. [53] Scalar-relativistic Kohn-Sham equations within the local 
density approximation [54] were solved taking all relativistic effects into account 
except for the spin-orbit coupling. The calculations were performed on 72 irreducible 
k points within the primitive wedge of the Brillouin zone. LMTO electronic structures 
were analyzed by calculating crystal orbital Hamiltonian populations (COHPs) and 
electron localization functions (ELFs). The crystal orbital Hamiltonian population 
COHP[55, 56] is a tool which permits energy-resolved analysis of bonding. The ELF 
provides a measure of the local influence of the Pauli repulsion on the behavior of 
electrons and permits the mapping in real space of core, bonding, and non-bonding 
regions in a crystal. [49, 50] Calculations were performed on the ideal crystal structure 
of pyrochlore Bi2Ti207[33] (space group Fd3m, a = 10.37949 A), displayed in figure 1. 
The atom positions are Bi (0,0,0), Ti (5,5,5), O (g,|,x) and O' (|,|,|). x = 0.43128 was 
taken from the reported room temperature crystal structure refinement. [33] LMTO 
calculations were also performed on the Aliam crystal structure of ferroelectric 
SrBi2Ta209 taken from the detailed structural study of Rae et al\b7\ A grid of 72 
irreducible k points within the primitive wedge of the Brillouin zone was used 
for the calculation on SrBi2Ta209. The crystal structure of the defect pyrochlore 
Pb2Sn206[48] is reported in the International Crystal Structure Database (entry 
15308) to be a cubic pyrochlore with a = 10.719 A. Atoms are Pb (0,0,0), Sn {\,\,\), 
and O {l,l,x) with x = 0.4250. O' is missing in this structure. LMTO calculations on 
Pb2Sn206 employed a grid of 72 irreducible k points. 
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Figure 2. (Color) (a) Densities of Bi s and p states for the ideal cubic structure of 
Bi2Ti207. (b) Densities of the p states of O and O' in Bi2Ti207. COHPs showing aU 
Bi-O and Bi-O' interactions within the primitive unit cell of Bi2Ti207. The origin on 
the energy axis is the top of the valence band. 



3. Results 

3.1. Bi2Tii07 

The Bi20' network of Bi2Ti207 can be imagined as being formed of a diamond lattice 
of O' , with one Bi atom at the center of every O'-O' linkage. The Bi20' network thus 
satisfies one of the Bernal-Fowler rules[ll] in that each O' is tetrahedrally surrounded 
by Bi. The network does not satisfy the "two-near and two-far" rule, meaning that 
ideal pyrochlore Bi2Ti207 is perfectly crystalline and has no residual entropy. 

The projected Bi and oxygen LMTO densities of state (DOS), and the Bi-oxygen 
COHPs of cubic Bi2Ti207 are displayed in the three panels of figure 2. The band gap 
in cubic Bi2Ti207 is about 2 eV. Bi s states in figure 2(a) are seen to reside in a very 
narrow window of energy around -10 eV with respect to the top of the valence band 
(set as the origin of the energy axis in this and other plots). Bi p states, as expected 
for a formally Bi'^+ system, are largely empty and form the conduction band which 
starts at 2 eV. Some Bi p states are found mixed with the valence band as a result of 
covalency with oxygen. Bi s and p states of Bi2Ti207 are resemble those previously 
observed in cubic "pre-distorted" BiMn03.[39] O p and O' p states are displayed in 
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Figure 3. (Color) Bi-O' network of Bi2Ti207 showing the valence electron localization 
isosurfaces visualized for a value of ELF = 0.65. The total ELF for a single plane is 
projected on the back of the unit cell. The localization scale runs from deep blue (ELF 
= 0) to white (ELF = 1). 

figure 2(b). There are six O for every one O', giving rise to a large difference in the 
relative number of states. O' states are particularly prominent around -5.5 eV, near 
the bottom of the valence band, where O' p forms band with Bi p states. Bi-O and Bi- 
O' COHPs are displayed per primitive unit cell (or for a single BiO(,02 polyhedron) 
in figure 2(c). The integrated COHP which is indicative of the strength of interaction, 
confirms that Bi s states are quite inert in the cubic crystal structure. Stronger Bi p 
and O' p interactions manifest as a large bonding COHP at the bottom of the valence 
band. Because of the longer distance between Bi and O (2.635 A) compared with Bi 
and O' (2.247 A), the Bi-O interaction is weak and O obtain most of their bond valence 
from Ti. 

The valence electron localization function (ELF) of Bi2Ti207 is plotted within the 
real space of the crystal structure in figure 3. For clarity, Ti and O atoms are not 
displayed. Bi and O' are connected to emphasize the "interrupted" diamond lattice. 
Isosurfaces of the ELF for a value of 0.65 manifest the lobes of s electrons around 
Bi. These lobes are constrained by the symmetry of the structure to form cylindrical 
objects with the cylindrical axis along different [111] directions, corresponding to the 
Bi-O' bonds. The ELF is also displayed on the back plane of the crystal structure as a 
map running from deep blue (poorly localized regions) to white (strongly localized 
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Figure 4. (Color) Crystal structure Aurivillius phase SrBi2Ta209 projected down the b 
axis. The black spheres are Bi, orange are O, grey are Sr and blue octahedra surround 
Ta. 

regions). The particular manner in which the lone pairs localize around Bi atoms 
is an artifact of the ideal crystal structure is used here. In reality, both Bi and 
O' are shifted from their crystallographic sites, [33] albeit in a incoherent manner. 
To understand what the lone pair localization should more properly resemble, the 
disposition of lone pairs in the Aurivillius phase ferroelectric SrBi2Ta209 is examined 
in the following subsection. 

3.2. SrBi2Ta20g 

The A2iam crystal structure of SrBi2Ta209 [57] is displayed in figure 4 and comprises 
alternate stacks of double-perovskite [SrTa207]^ slabs interleaved with [Bi202]^^ 
slabs. The latter are reminiscent of the Bi-O sublattice of BiOF[41] where Bi lone 
pairs organize themselves in an antiferrodistortive manner, akin to what is seen in 
PbO.[40, 41] In this distorted crystal structure, Bi s states [figure 5(a)] are significantly 
more disperse than was observed for cubic Bi2Ti207. Relatively few Bi p states are 
observed in the region of filled O p states [figure 5(b)]; most Bi p states are empty 
and, starting at 2 eV, form the conduction band along with empty Ta d states (not 
displayed). The Bi-O COHP displays quite strong bonding in the region of Bi s states. 
Both Bi s states and O p states are filled, so their interaction gives rise to filled bonding 
states around -9 eV and filled antibonding states just below the top of the valence 
band. The intermediation of empty Bi p states provides net stabilization. [42] The Bi- 
O COHP of SrBi2Ta209 displayed in figure 5(c) is what the Bi-O COHP of Bi2Ti207 
displayed in figure 2(a) would resemble were it permitted to take on an appropriately 
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Figure 5. (Color) (a) Densities of Bi s and p states in SrBi2Ta209. (b) Densities of the all 
p states of O SrBi2Ta209. COHPs showing all Bi-O interactions within the primitive 
unit cell of SrBi2Ta209. The origin on the energy axis is the top of the valence band. 



distorted coordination. 

The valence ELF of SrBi2Ta209 is displayed in the different panels of figure 6. 
The left panel of figure 6 displays the ELF isosurf ace within the space of the unit cell. 
For clarity, other atoms than in the Bi-O network have been let out of the depiction. 
The net "polarization" of the lone pair lobes is more clearly seen in projections down 
the c axis [figure 6(a-d)]. The ELF was been calculated on a 48x48x144 (a x b x c) 
real-space grid within the unit cell. Figure 6(a-d) shows four sections of the ELF at 
four different heights along the c direction, corresponding to positions of Bi atoms 
in the unit cell. In this view, it is clear that all the Bi lone pairs have a component 
that points in the a direction, meaning that the Bi are shifted slightly along b. It is 
this in-plane shift which gives SrBi2Ta209 its polarization. The very different shapes 
of the lone pair lobes associated with centered Bi in Bi2Ti207 (figure 3) and with off- 
centered Bi in SrBi2Ta209 (figure 6) should be noted. 



3.3. Pb2Sn20e 

Morgenstern-Badarau and Michel[48] have reported the defect pyrochlore crystal 
structure of "Pb2Sn206 ■ x(H20)". Interestingly, the compound has long history and 
could be one of the components of a yellow pigment used by renaissance painters 




Figure 6. (Color) ELF isosurface projected within the unit cell of SrBi2Ta209 for an 
ELF value of 0.65. The panels (a-d) on the right are planes of the ELF parallel to the 
c axis; the planes are chosen to contain Bi atoms. The values of z/ c of the different 
planes are indicated as ratios. 



called lead-tin yellow. [58] The O' atom is missing in this crystal structure and is 
replaced by a vacancy (□). A2O' in the usual pyrochlore then corresponds to Pb2n. It 
is natural to assume that Pb^+ lone pairs would occupy the site that is normally taken 
up by O'. This would parallel the crystal-chemical relation between litharge PbO and 
Pb02; the former having lone pairs in the place of the oxygen atoms of the latter. [59] 
However the A2O' network of the pyrochlore lattice depicted in figure 1 allows us to 
recognize that for every four Pb atoms, there are two sites into which the lone pairs 
can localize. The problem now maps precisely on to the Bemal-Fowler rules with 
the lone pair-Pb connection corresponding to the short (covalent) O-H interaction. 
The valence ELF of cubic Pb2Sn206 displayed in figure 7 support this argument. 
The symmetry of the crystal structure prevents localization into a single lobe, and 
instead, a bi-lobed lone pair is obtained. Topological frustration again prevents 
coherent ordering of the lone pairs in some distorted low-temperature structure. 
Re-examination of the crystal structure of defect pyrochlore Pb2Sn206 should reveal 
disordering on this site; splitting every Pb atom into two along the D-Pb-D direction 
(the [111] directions). The precise nature of such splitting would therefore be distinct 
from what is observed for Bi in Bi2Ti207 where it is forms an annulus around O'-Bi- 
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Figure 7. (Color) ELF isosurface projected within the unit cell of Pb2Sn206 for an ELF 
value of 0.65. The yellow spheres are vacant positions corresponding to O' in the ideal 
pyrochlore. 

0'.[33] 

4. Implications 

To aid the rest of the discussion, the panels of figure 8 display two different views of 
the O'-centered tetrahedral network, also referred to as the tridimyte sublattice,[6] of 
the pyrochlore structure. The two views suggest that we could view the pyrochlore 
structure as a stellation followed by stacking of kagome lattices of A atoms. A kagome 
lattice is a lattice of corner sharing triangles. The stellation of this lattice is achieved 
by making each triangle a tetrahedron by alternately capping the triangles above 
and below the planes they are in and creating tetrahedra pointing either up or down 
(along [111] in the pyrochlore structure). In 2D the pyrochlore lattice should be well- 
represented by a kagome lattice. This allows us to compare in a simple manner, the 
topologies of perovskites, which we represent by a square lattice, and the pyrochlore, 
which we represent using the kagome lattice (figure 9). 

A square lattice is easily distorted to a rectangular lattice [figure 9(a) and 
(b)] allowing coherent localization of lone pair lobes, which we indicate with 
arrows. This simple description is reasonably accurate in describing the pyroelectric 
to ferroelectric (cubic Pm3m to tetragonal P4mm) phase transition in perovskite 
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Figure 8. (Color) (a) Corner-connected 0'A4 tetrahedra of the A2B2O6O' structure 
projected down [111] (b) A single slab of the network seen from the side; the 
pyrochlore structure described as a stellated kagome network. 





Figure 9. (Color) (a) is a square lattice with lone-pair active ions in the corners. 
Coherent localization of the lone pairs results in a polar rectangular lattice (b). (c) 
is a kagome lattice with lone pair active ions at the corners. The lone pairs can either 
localize towards other vertices (d) or towards the centers of the triangle (e), but they 
do not do so coherently. 
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PbTi03.[18, 41] 

A kagome lattice with lone pair-active ions at its corners [figure 9(a)] could 
localize lone pairs in two different ways: The first is for the lone pairs to localize 
perpendicular to the line joining the centers of the triangles, as depicted in figure 9(c). 
This is the case with the annular lone-pair localization in Bi2Ti207. The second is for 
the lone pairs to localize towards the centers of triangles, as depicted in figure 9(d). 
This is representative of kagome spin ice[60] and is similar to what is observed 
Pb2Sn206. In neither case is it a simple matter to arrange the arrows such that a unit 
cell is obtained. Note that neither of these arrangements would he ferrodistortive. In 
fact, while we have focused on the paucity of ferroelectric pyrochlores, our analysis 
actually refers to any symmetry lowering phase transition. 

In the rare instances when pyrochlore oxides distort as a result of coherent 
ordering of stereochemically active lone pairs, the resulting crystal structures 
can be be quite complex, for example, as determined recently for Bi2Sn207.[61] 
However, retention of the cubic structure is important for properties. For 
instance, the pyrochlore composition Bi2Zni/3Nb4/307 crystallizes in a monoclinic 
structure related to zirconolite,[62] and the material is not as interesting as 
(Bii 5Zno.5)(Nbi.5Zno.5)07 for high-fc applications. Wang et al.[25] have presented a 
simple and powerful argument for when A2B2O7 pyrochlores with lone pair active 
A cations distort. They point to a relative competition between covalency in the A-O 
and B-O networks. Pyrochlores with strongly covalent B-O networks tend to retain 
the cubic structures. 

The arguments presented here bear close similarities to the frustration of soft 
modes in ZrW208, proposed by Ramirez and coworkers[65] as being responsible 
for its negative thermal expansion behavior. [66] Recently, Loidl and coworkers[67] 
reported unusual magnetoferroelectric properties of spinel CdCr2S4. In this 
compound, relaxor ferroelectricity is observed in the absence of any substitutional 
disorder. The authors suggest that the origin of such relaxor behavior could be the 
topological frustration of a true ferroelectric ground state. In the superconducting 
/3-pyrochlore oxides, Kunes, Jeong, and Pickett[68] have used first principles 
calculations to determine that local displacements of the A cation are critical in 
determining superconducting TcS. One could argue that in the absence of the 
topological frustration provided by the pyrochlore lattice, these local displacements 
would lead to a distorted ground state structure. In a related vein, for macroscopic 
materials built up from bars and pins, the unusual properties of kagome lattices with 
respect to the applications of stresses have already been noted. [69] 

To summarize, it is proposed that in pyrochlores such as (Bii,5Zno.5)(Nbi.5Zno.5)07, 
properties are dominated by topological frustration rather than cation site disorder. 
Bi2Ti207 itself is therefore a candidate for high-fc behavior, and should be regarded 
as a parallel of systems such as cubic, predistorted PbTiOs. In relaxor ferroelectrics, 
substitutional disorder frustrates long-range ordering and gives rise to unusual be- 
havior, such as broad dielectric anomalies as a function of temperature. [63, 64] If re- 
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laxor ferroelectrics with site disorder are the ferroelectric equivalent of substitutional 
spin glasses such as Cu:Mn, then pyrochlore dielectrics could be considered the fer- 
roelectric analogues of topologically frustrated magnets. In analogy with topologi- 
cally frustrated magnets, the recording of low-temperature heat capacity of Bi2Ti207 
should permit the residual entropy to be measured and compared with values found 
in pyrochlore spin-ice. [15] Finally, in analogy with spin-ice, it is proposed that the 
term "charge-ice" is used to describe frustrated pyroelectric phenomena in systems 
such as spinels and pyrochlores. Topological frustration in polar materials could con- 
stitute a powerful design parameter for useful and novel materials properties. 
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